Quasistatic Stick-slip in Dislocation Core and Prenkel-Kontorova Chain 
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By means of atomistic simulations, we demonstrate that a dislocation core exhibits intermittent 
quasistatic restructuring during incremental shear within the same Peierls valley This can be re- 
garded as a stick-slip transition, which is also reproduced for a one-dimensional Frenkel-Kontorova 
chain under rigid boundary conditions. This occurs due to a discontinuous jump in an order pa- 
rameter of the system, which signifies the extent of region forbidden for the presence of particles in 
the chain. The stick-slip phenomenon observed in the dislocation core is also shown to be reflected 
after dimensionality reduction of the multidimensional atomic coordinates, which provides a basis 
for comparison with the simple one-dimensional chain. 

PACS numbers: 61.72.Lk, 02.70.Ns, 45.05. +x 



How much illustrative can we expect an ideal system 
to be, which is regarded as the simplified model of a real 
physical phenomenon? The answer to this question is rel- 
evant in view of the successful applications of such models 
throughout the history of physics. The Kronig Penney 
model of band structure, Ising model for magnetic sys- 
tems, liquid drop model of atomic nucleus, cellular au- 
tomata model of self-organized criticality etc. are some of 
the revolutionary models with remarkable achievements. 
The Frenkel Kontorova (FK) chain [l] [2] is one such cele- 
brated model dealing with the mechanics of discrete non- 
linear systems. Conventionally, this model has success- 
fully been employed in studying a wide range of physical 
phenomena like the physics of dislocations [3] and crow- 
dions [4 in metals, adsorption of atoms on crystal sur- 
face [5 , magnetic structures [6] etc. In addition, it has 
also rendered rich understanding of processes like col- 
loidal friction [7] and biopolymers [3 [9] in recent years. 

The FK model consists of a linear chain of particles 
connected by springs and placed over a substrate poten- 
tial. In its most elementary form, the springs are assumed 
to be Hookean, while a sinusoidal substrate potential is 
considered. It was first envisaged [lj to represent the 
dislocation core as a kink in the FK chain. In particu- 
lar, the existence of a threshold force to move the kink 
directly corresponds to the Peierls stress [10] for the dis- 
location motion. The Peierls stress for a dislocation is 
often computed by means of atomistic simulation [TT1 - 
H5] , where the incremental shear stress (or strain) is ap- 
plied quasistatically to the crystal at T = K until the 
dislocation moves to the next lattice site from the previ- 
ous one. Similar quasistatic simulation of the FK chain 
[2j [14] also yields the minimum force required to move 
the chain. Thus, the FK model thematically depicts the 
mechanism of crossing the Peierls barrier in perspective 
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of discreteness and nonlinearity intrinsic to the lattice. 
Nevertheless, it is obscure whether this simplified model 
is capable of providing finer details of the process. A 
closer look into this issue is worthwhile as it would not 
only highlight the extent of resemblance between the re- 
alistic physical system and the representative model, but 
also enable us to extend the applicability of the model 
beyond a coarse representation. In this Letter, we study 
the atomistic simulation of forcing a dislocation core out 
of its Peierls valley to obtain the atomic trajectories at 
sub-Burgers vector resolution. Interestingly, the simu- 
lations reveal the occurrence of intermittent relaxation 
bursts at such fine scale. Aiming to unfold the underly- 
ing mechanism, we apply the appropriate boundary con- 
ditions to a simple one-dimensional FK chain to let it 
mimic the boundary conditions used in the atomistic sim- 
ulation. Even though the FK model is regarded merely 
as a conceptual tool subsuming the inherent discreteness 
and nonlinearity of the lattice, it also exhibits the abrupt 
bursts of structural relaxation, similar to the dislocation 
core. In what follows, we shall argue that such a striking 
feature can be perceived as the quasistatic counterpart 
of the stick-slip motion frequently observed in numerous 
dynamic systems [T5HT8] . In addition, the technique of 
principal component analysis (PCA) [19] has been used 
in an innovative way to establish the correspondence be- 
tween a real physical system and its ideal model of lower 
dimensionality. 

The Peierls stresses, which are defined at absolute zero 
temperature, have been computed here [20] for the four 
metals, molybdenum, iron, aluminium and copper. The 
simulation scheme is akin to that used in the earlier mea- 
surements [TTJ [13] . An edge dislocation is introduced in 
a slab of finite thickness with periodic boundaries along 
the directions of dislocation line and Burgers vector. The 
x, y and z dimensions of the b.c.c. simulation cells are 
90.5a(lll)/2 , 40a(101) and 5a(121), where a is the lat- 
tice constant. The Burgers vectors and the dislocation 
lines are along the x and z directions respectively [see the 
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FIG. 1: (color online), (a) Rise in potential energy of the 
Mo crystal with a typical simulation cell in the inset. Note 
the sudden drop marking the instant of crossover to the ad- 
jacent Peierls valley. Differential displacement profiles for the 
dislocation core atoms in (b) Mo, (c) Fe, (d) Al and (e) Cu. 



inset of Fig. 1(a)]. Corresponding cell dimensions for the 
f.c.c. systems are 90.5a(101)/2 x 20a(lll) x 5a(121), 
where the perfect dislocations splits into Schockley par- 
tials on account of the Frank's criterion [25 . Shear strain 
in the system is increased in small steps by tilting the ver- 
tical boundary of the simulation cell. At each step, the 
system is relaxed to the minimum energy configuration 
using the conjugate gradient method [11], while the top 
and bottom surfaces are kept fixed during the relaxation. 
The potential energy of the system is recorded and the 
crossover of the dislocation core to the next Peierls val- 
ley is marked by a drastic drop in the total potential 
energy profile. Detailed analysis of the simulation out- 
put involves the region only up to the crossover points, as 
shown in Fig. 1(a). The dislocation core atoms are iden- 
tified by using suitable centrosymmetric deviation pa- 
rameter (CSD) windows [11]. For b.c.c Mo and Fe, the 
CSD window of range 1.4-10 A 2 has been used, while the 
ranges 3-20 A 2 and 3-16 A 2 are employed for Al and Cu 
respectively [2T] . 

As the process of applying incremental shear strain 
is followed by relaxation, the atomic structure tends to 
reconfigure so that the potential energy of the system is 
minimized. The core structure can be specified by a set of 
n c vectors, say {r^} (i = 1,2,.., n c ), where N c is the num- 
ber of atoms in the dislocation core. We can now quantify 



the extent of aggregate core displacement at the n ttl step 
of incremental shear strain with respect to the previous 

step as A = yX^=i \ r i( n ) ~ r i( n ~ 1)| 2 - These differ- 
ential displacements of the dislocation cores are plotted 
with the shear strains in Figs, l(b-e). Surprisingly, one 
can identify the intermittent relaxation bursts character- 
ized by narrow peaks in the profiles. Obviously, a small 
value of A means that r^(n) « r^(n — 1) implying only 
a small structural change in the dislocation core from 
the previous step of applied strain, whereas a large value 
is indicative of drastic structural rearrangement. This 
proves that instead of exhibiting a continuous response 
to the incremental strain, the dislocation core structure 
remains almost locked in between two successive bursts. 
This quasistatic phenomenon is apparently analogous to 
the stick and slip states of dynamic variables in many 
physical processes of interest. 

In this study, the b.c.c. metals Mo and Fe are simu- 
lated using the modified Finnis-Sinclair interatomic po- 
tential [22], whereas the glue potential [23] and an em- 
bedded atom model [24] are used for Al and Cu respec- 
tively. In spite of a variety of potentials used in these 
computations, the relaxation bursts can be observed in 
all the studied systems. This suggests that the ori- 
gin of these bursts lies rather in a more fundamental 
physical phenomenon, than in the complexities of the 
potential models employed in the simulations. In this 
context, the possibility of closely investigating the FK 
model arises as it reflects the basic features of a dis- 
location core, namely, the nonlinear nature of interac- 
tions and discreteness of the lattice. The potential en- 
ergy of a finite FK chain with N connected particles is 



u = Er=i [f(*i+i 



l) 2 + E(l-cos*f±)], where 



the first term denotes the harmonic spring potential with 
Xi, the position of i th particle, /, the equilibrium length 
of each spring and ft, the spring constant. The second 
term represents the periodic substrate potential where b 
and E are the periodicity and magnitude of the poten- 
tial respectively. Here we assume the natural length of 
the springs connecting the adjacent atoms to be equal to 
the periodicity of the substrate potential (i.e., I = b = 
1 arb. unit). At this point, it should be noted that de- 
spite its success in demonstrating the origin of the Peierls 
barrier, the FK chain does not exhibit a direct correspon- 
dence to the dislocation core in all aspects. For instance, 
the long range interactions between two dislocations [25] 
is in clear contrast to the exponential short range kink- 
kink interactions in the FK chain [2]. Thus, a judicious 
choice of boundary conditions and parameters of the FK 
chain is necessary so that the system can distinctly show 
the features of interest. In the atomistic simulations, 
the top and bottom surfaces of the crystalline slabs were 
kept fixed during the relaxation process so that the sys- 
tem could not revert back to the previous state of shear 
strain. To implement this on our FK chain, we fix both 
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FIG. 2: (color online), (a) Differential displacement and (b) 
change in the length of springs for the FK chain accommodat- 
ing N — 91 particles in N v = 60 valleys with E/k — 0.01. (c) 
At E/k = 0.035, the differential displacement profile shows in- 
termittent peaks (P1-P8) along with three small humps (HI, 
H2, H3). (d-g) Changes in the lengths of the 90 springs show- 
ing the transition from the stick states (d,e) to the slip state 
(f) at the fourth relaxation burst (P4), followed by another 
stick state (g) corresponding to the instants indicated by the 
four arrows in (c). (h) Variation in the order parameter ^ as 
a function of shift in the chain, $. The features corresponding 
to the humps and peaks in (c) are marked here as well. The 
peaks in (c) always coincide with discontinuous jumps in ^. 



of its ends, thereby yielding a fixed length and fixed den- 
sity condition. We accommodate N particles in N v val- 
leys so that the coverage parameter N/N v is close to ~ | 
to approximately resemble the local coverage due to the 
extra half-plane of atoms at the core of an edge dislo- 
cation. The chain is gradually shifted in small steps of 
S and relaxed after each shift keeping the first and the 
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FIG. 3: (color online). Schematic representation of discontin- 
uous rise in the order parameter ^ as the sudden broadening 
of the forbidden region. This causes the particle in the FK 
chain to abruptly cross the potential hill. 



last particles rigid during the relaxation process. Hence, 
the net shift of the chain at the n th step can be rep- 
resented by the coordinates of the fixed end particle as 
& = x\ = n5, where the incremental shift (S) should 
have been infinitesimally small for the ideal quasistatic 
process. However, due to the trade-off between compu- 
tational time and resolution, a small value of 5 = 10 -2 
(arb. units) is chosen and found to be sufficient to pro- 
duce the requisite spatial resolution. Moreover, in the 
present case of smooth periodic substrate potential, the 
simple steepest-descent algorithm [11] reasonably yields 
the relaxed states. 

The differential displacements are now computed for 
the particles of the FK chain for different E/k ranging 
from 0.01 to 0.1. Fig. 2(a) shows the differential dis- 
placement for E/k = 0.01 as a function of the shift, <£, 
imparted to the chain where a continuous wavy nature 
can be observed. In addition, the change in the lengths 
of the connecting springs in between two successive steps 
is shown in Fig. 2(b), where the change is noticed along 
the entire chain. With increase in the value of E/k, this 
wavy nature changes to relaxation bursts characterized 
by the intermittent peaks, which can be seen from Fig. 
2(c) for E/k = 0.035. Remarkably, these peaks which 
demonstrate the phenomenon of quasistatic stick-slip are 
present here also, similar to those in Figs, l(b-e). Cor- 
responding changes in the length of springs, which are 
confined to the regions near the fixed ends of the chain 
[Fig. 2(f)] are also intermittent and synchronized with 
the occurrence of relaxation bursts. 

For the FK chain in its ground state, the minimum dis- 
tance of the particle from the nearest maximum of the 
substrate potential has been used as an order parame- 
ter by Coppersmith and Fisher [26] to characterize the 
Aubry-transition and breaking of analyticity [2 . For the 
plot shown in Fig. 2(c), E/k is large enough to cause the 
breaking of analyticity. However, in this study we en- 
counter another variable <£>, denoting the extent of shift 
in the chain, which determines the ground state config- 
uration of the entire chain. Therefore, the order param- 
eter given by \I/ = min\xi(mod 1) — 0.5|^i 5 jv has been 
computed and displayed in Fig. 2(h) as a function of 
<£. A finite non-zero value of \I/ signifies the span of re- 



gion around the top of the substrate potential where the 
presence of a particle is forbidden. It is noticeable that 
despite the apparently sudden occurrences of relaxation 
bursts, the instant of peak in Fig. 2(c) is always preceded 
by a gradual drop in \I/ [Fig. 2(h)] during the stick state, 
thereby indicating the gradual narrowing of the forbidden 
region. Thereafter a discontinuous jump in ^ coincides 
with the transition to slip state at which the relaxation 
burst occurs. This is indicative of abrupt broadening 
of the forbidden regions and consequently, we expect at 
least one particle to suddenly cross over the peak of the 
substrate potential (see the schematic in Fig. 3) thereby 
causing a steep rise in the differential displacement pro- 
file. In addition, Fig. 2(h) also shows three instances 
where \I/ drops and reaches values close to zero, and then 
rises continuously. Because of this continuous change in 
the order parameter, we observe small humps in the dif- 
ferential displacement profile which is in sharp contrast 
to the abrupt peaks, where the rise in ^ is discontinuous. 
This can be somewhat difficult to identify, for example, 
during the second hump (H2), where it coincides with 
occurrence of fifth peak (P5) in the observed profile. 

Typical ground state configurations of the FK chain is 
given in Fig. 4(a), where most of the particles follow a 
structural pattern associated with a highly stable energy 
state and attempt to maintain it despite the incremen- 
tal change in <£. As a result, the local configurations 
near the two fixed ends is out of skew with the rest of 
the chain and total energy of the system increases un- 
til the next relaxation burst. During this stick state the 
differential displacement of particles are found to be neg- 
ligibly small and hence, the coordinates of the particles, 
Xi(mod 1) after relaxation, when plotted with respect 
to the coordinates before energy minimization show the 
linear behavior [Fig. 4(b)]. However, during a relax- 
ation burst some of the coordinates are knocked out of 
the straight line pattern, as shown in Fig. 4(c). As pre- 
dicted earlier, there is always one or more particles which 
cross the peak of the substrate potential corresponding 
to Xi(mod 1) = 0.5 as marked in the figure. 

The FK chain presented here is a 1-D model of an es- 
sentially 3-D atomic structure of a dislocation in crystal. 
Each and every relaxation burst always drives the chain 
in a forward direction as both ends are shifted in same 
direction. Similarly, it must be ascertained that the re- 
laxation bursts as observed in the atomistic simulation of 
dislocation core also drive the core in an effectively for- 
ward direction. The notion of directionality is trivial in 
the FK model because of the inherent single dimension- 
ality associated with the structure. However, the motion 
of a dislocation can be perceived only on a coarse scale 
of length, where the dislocation hops from one lattice 
site to another. On a sub-Burgers vector scale of length, 
it is difficult to associate a sense of unidirectional mo- 
tion to the self-assembly of core atoms within the same 
Peierls valley unless a specific directionality is ascribed 
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FIG. 4: (color online), (a) Two typical ground state config- 
urations of the FK chain ($=0, 0.9) with N = 91 particles 
in N v = 60 valleys for E/k = 0.035. The relaxed coordinates 
of the particles (b) just before and (c) at the instant of oc- 
currence of the peak P5 in Fig. 2(c) plotted with respect to 
the unrelaxed coordinates. The particles which have crossed 
the maxima of the substrate potential i.e., Xi(mod 1) = 0.5 
(shown by dotted horizontal line), are identified in the rect- 
angular frame. 



to it. Clearly, the conventional technique of describing 
the core position as the center of mass of all the core 
atoms (Refs. [TTJ [27j [28] for example) lacks the requisite 
resolution and alternative data mining tool needs to be 
explored. In this scenario, we opt to use the PCA [19] 
as a prolific tool capable of providing a high degree of 
compressibility of a high dimensional data and furnish 
its projection on a hyperspace of reduced effective di- 
mensionality. The versatility of the PCA is reflected in 
its successful applications across a wide range of studies 
[29H33] etc. In the present case of atomistic simulations, 
the coordinates of core atoms are recorded at each step 
of incremental shear strain and this strain series data is 
arranged as a n s x 3n c matrix, where n s denotes the num- 
ber of strain steps and n c is the number of core atoms. 
Now each of the 3n c columns is separately mean cen- 
tered [34 and the mean-deviation matrix thus formed is 
used to generate the covariance matrix [19]. Diagonal- 
ization of the covariance matrix yields the eigenvalues, 
and the corresponding eigenvectors. Interestingly, for all 
the metals under study, the largest normalized eigenval- 
ues are always found to be in excess of 90%. Such large 
values conclusively prove a high compressibility intrinsic 
to these sets of multi-dimensional data. The projections 
of the datasets along the principal directions correspond- 
ing to the largest eigenvalues are presented in Fig. 5 for 
Mo and Al, for example. The sudden jumps present in 
these devil's staircase-like profiles are synchronized with 
the occurrence of peaks in Fig. 1(b) and (d) and are typ- 
ical signatures of the stick-slip process (see the first and 
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FIG. 5: (color online). Principal projections for (a) molyb- 
denum and (b) aluminium show staircase-like profiles. The 
quasiplateaus represent the stick states whereas the slip states 
are reflected as the sudden jumps. 



second panels of Fig. 4 in Ref. [9 ). Moreover, one can 
also observe that each monotonic jump in the projected 
profile always causes a translation in the same direction, 
thereby offering a ground for comparison with the FK 
model. 

To conclude, we have shown that at sub-Burgers 
vector resolution, intermittent relaxation bursts occur 
in the quasistatic simulation of dislocation core. Similar 
features are observed for the simple one-dimentional FK 
chain as well. This is attributed to a transition of the 
system from an effective stick to slip state on account of 
abrupt broadening of the forbidden region around the 
peak of the substrate potential. During the stick-slip 
transition, one or more atoms cross maximum of the 
substrate potential and hops over the forbidden zone to 
cause a prominent rise in the differential displacement. 
Moreover, the tool of principle component analysis has 
been used in an innovative way to extract the effective 
dimensionality of the atomistic data of the dislocation 
core atoms. The projections of the atomic trajectories 
on the principle directions further corroborate the 
efficacy of the 1-D FK chain in revealing the complex 
3-D structure of the dislocation core. 

A. Dutta acknowledges the financial support from 
CSIR, India to carry out this research work. 
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